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Abstract: We calculate the running of the quark mass and the quark condensate using a djmamical 
soft-wall model by Csaki and Reece. We show that a correct running can be obtained with an appro- 
priate sign for the dilaton field. In the soft-wall model with quadratic dilaton profile, a wrong sign 
for the dilaton field can give rise to a massless vector meson, a signature of spontaneously broken 
symmetry in the vector sector, which is not observed in nature. With a right parameter range, we ob- 
tain no such massless vector meson. We also find that, contrary to the soft-wall model with quadratic 
dilaton profile, the model allows a non-vanishing quark condensate in the chiral limit. We obtain the 
vector meson mass spectrum similar to that of the hard-wall model. By varying a parameter in the 
model we can fit the running of the quark mass and of the quark condensate. 
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1 Introduction 

Understanding Quantum Chromod5mamics (QCD) at the low-energy region has been a challenge for 
physicists in the past half century. While the high-energy region of QCD is asymptotically free [1] 
and can be studied and understood well using perturbative method, in the low-energy region, QCD 
interactions actually become strong and non-linear and cannot be studied perturbatively. Another 
challenge is to understand the structure of nucleons as the non-linearity of QCD cause the interactions 
between quarks in the low-energy region to be very complicated. 

For the lack of formal derivation, most of the efforts to understand QCD at low-energy region 
are concentrated in building models that closely resemble QCD at this energy scale, in particular in 
modeling hadrons. Many models have been proposed in the past, such as Bag Model [2], Constituent 
Quark Model [3], Chiral Soliton Model [4], Chiral Lagrangian [5, 6], etc. While most of these models 
work quite well in the energy region they are prescribed and explain certain properties and behavior 
of hadrons, they are not derived from first principle. They also generally cannot extend to the high- 
energy region where QCD works very well. 

In 1997 Maldacena conjectured that a certain type of string theory, in particular the type IIB string 
theory in AdSs x S5 corresponds to the conf ormally invariant large Nc SU(Nc) A/" = 4 Supersymmetric 
theory on the 4-dimensional space on the surface of the 5-dimensional AdS space [7]. This is how the 
moniker AdS/CFT got its name. With this conjecture, the relationship becomes more tangible and 
a dictionary can be composed to relate parameters and symmetries in one theory to another [8, 9]. 
Better yet, the two regions of the theory cross where strong coupling, i.e. non-linear region in one 
theory corresponds to weak coupling, i.e. perturbative region in another theory. To calculate a theory 
in the strong region all we have to do is to calculate the corresponding theory in the perturbative 
region. There are two approaches in constructing the dual of QCD, dubbed the top-down and the 
bottom-up approaches. In the former, one starts from a string theory with an appropriate brane 
configuration [10-12]. In the latter, one begins by introducing a 5D action with the few fields needed 
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to reproduce the well-known properties of low-lying QCD states [13-15]. In this paper we will follow 
the bottom-up approach. 

In its simplest form, the bottom-up models use a 5D AdS metric with cutoff to break conformal 
symmetry and simulate confinement in QCD. The 5D fields present in the action are bifundamental 
scalar field X{x,z) and gauge fields, L(x, z) and R{x,z), satisfying the SU(2)lX SU(2)r symmetry. 
These 5D fields corresponds to qiq^ix), Jl{x) and }r{x) operators in the 4D gauge theory respec- 
tively. The model produces QCD features quite well [16-24]. However, the model has some draw- 
backs. There are ambiguities in the infrared boundary conditions, and the resulting meson spectrum 
does not satisfy the observed Regge trajectory where the mass squared grows as n. Instead, the spec- 
trum of the mass is similar to that of the quantum system of a particle in a box, i.e. grows as 

Another way to simulate confinement is by introducing a background dilaton field to the 5D 
action [25]. In order to obtain a meson spectrum that satisfies linear Regge trajectory, the profile of 
the dilaton field is chosen to be quadratic while the metric of the 5D space remains AdS. However, 
while the soft-wall model successfully reproduce the linear Regge trajectory, it fails to reproduce a 
QCD-like chiral symmetry breaking. In particular, in the limit of vanishing quark mass parameter, 
the quark condensate parameter is vanishing as well. A way to overcome this problem, among other 
things, was pointed out in [26]. They introduced a quartic term to the bulk scalar potential, included 
a tachyon field and modified the dilaton profile. 

A dilaton field with opposite sign, relative to the one used in [25], has been considered in the 
literature [27-29] as it seems to have more appealing confinement properties. In [25] and later reem- 
phasized in [30], this reversed sign choice was discarded because it leads to an unphysical massless 
vector meson. It was shown that the problem can not be easily remedied imless the asymptotic value 
of the dilaton field at the IR boimdary is changed. With the correct dilaton sign, the wave ftinction of 
the massless mode is non-normalizable, hence, decouple. 

Various aspect of QCD has been investigated using the holographic principle [31-41]. In this pa- 
per, we will use a model introduced in [42] where back reaction by the dilaton field to the AdS metric 
is considered. Even in the simplest case in which no potential term for the dilaton field presents, 
the back reaction from the dilaton field automatically provides an IR cutoff. With an addition of ap- 
propriate potential terms for the dilaton field, the model can reproduce asymptotic freedom of QCD 
quite well. 

The organization of this paper is as follows. In section 2, we will review the original AdS /QCD 
proposal for the hard-wall and the soft-wall model. In section 3, we will review the holographic 
model in [42]. We will discuss the vector part of the action in section 4, and show that the spectrum 
does not contain massless vector meson. In section 5, we consider the scalar part of the action and 
investigate the running of the quark mass and the quark condensate. Finally in section 6, we provide 
the conclusions. 



In AdS/CFT there is a correspondence between 5-dimensional fields (p{x,z), and 4-dimensional op- 
erators 0{x) living on the boundary of the 5-dimensional Anti-de Sitter (AdSs) space, where z is the 
fifth coordinate. The metric of AdSs can be written as 



2 AdS/QCD 




(2.1) 
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The fifth coordinate, 2, corresponds to energy scale (or momentum transfer Q^), i. e., z — > corre- 
sponds to infinite energy (Q^ — > 00). 

The 4D sources in the 4D generating functional Z4i5[^'^] is called (p^{x) and the generating func- 
tional is defined as follow: 

Zio[(p°] = ^exp (^ij . (2.2) 

The correspondence maybe written as 



-4D 



e'SsoIfe] ^ (2.3) 



where on the right hand side S^^lcpci] is the action evaluated for the classical solution (p^i{x,z) to the 
field equations with boundary condition 

lini</),,(x,z)=zV(^)- (2-4) 

The constant A depends on the nature of the operator O. Utilizing the correspondence, one obtain 
the n-point correlator 



(0|rO(xi)...O(x„)|0)= ^ 



3(xa).. 



(2.5) 



There are an infinite number of 5D fields corresponding to an infinite number of 4D opera- 
tors in QCD. However, there are only several crucial operators involved in the dynamics of spon- 
taneous chiral symmetry breaking. One of these crucial operators is the quark bilinear ipi^R, whose 
non-vanishrng vacuum expectation value breaks the full symmetry S\J {N j^) ixSU {N r down to 
S\J{Njr)v- The corresponding 5D field of the quark bilinear operator is a scalar field, X, which trans- 
form as a bifundamental. The conserved current operators of the SU(JV^)lxSU(N^)r symmetry, 
}1^{x) and Jr'{x), correspond to the 5D gauge field L^(x,z) and RI{x,z), respectively. 

2.1 Hard-wall 

In the hard- wall model [13], where the 5-dimensional space is AdS with cutoff, the meson part of the 
action can be written as 

d'x^r{\DX\+3\X\^ - l^{Fl + Fl)}, (2.6) 

with D^X = d^X - fAf X + iXA^ and F^^ = d^A^,^ - d^Af^j^ - i[Af^j^, A^,^]. The field X can 
be written as a product of a background scalar field Xo(z) and a pseudo-scalar field n''{x,z) 

X(x,z) = Xo(z) exp(/27r''f''). (2.7) 

In the case of Njr = 2 and assuming isospin symmetry, one can take Xg as multiple of the identity, and 
one does not have to worry about the ordering. The solution to the backgroimd scalar field equation 
of motion can be written as 

Xo{z) = ^(m^^z + ^CT^z^^ , (2.8) 

where the constant m,j can be identified as the quark mass sourcing the bilinear operator ipq^ij and 
the constant CTq as the vacuum expectation value of the quark bilinear operator. The constant ^ = 
y/WcIln was discussed in [43, 44]. 
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2.2 Soft-wall 



In the soft-wall model [25], confinement is modeled with a background dilaton field ^(z) = —kz^. 
This dilaton profile is chosen so that the meson spectrum behaves like, ~ n + S, for excitation 
number n and spin S. The dilaton field modifies the Lagrangian (2.6) in the form of an additional 
exponential factor e't'. The meson part of the action can be written as 

S= [d^xe^^r{\DX\+3\X\^-^{Fl + Fl)}, (2.9) 

The scalar field satisfies 

3^X0 - (iKz + d,Xo + ^Xo = 0, (2.10) 

with the following solution 

Xo(z) = CIZU Q,0,KZ2^ +C2Z^M (^|2,Kz2^ , (2.11) 

where U dan M are Kummer's functions. However, in order that the action evaluated on the solution 
is finite, the second term in the above solution must be dropped. Expanding the remaining term in 
the small z limit, one obtains 

Xo(z) = ^ z + ^ (l + 7£ + ln(KzV4)) z3 , (2.12) 

where 7e is the Euler-Mascheroni constant. Identifying c\ = mq^^/rz/'i, the above solution leads to 
a (X niq. In the limit niq 0, we eliminate both the explicit and the spontaneous symmetry breaking, 
in contradiction with QCD. 

In this paper, we use a model in which the dilaton field is dynamic. In the next section, follow- 
ing [42], we will choose a profile for the dilaton field that produces asymptotic freedom of QCD and 
obtain the potential term associated with the dilaton profile. We find that the model allows for a 
non- vanishing quark condensate in the limit of — > 0. 



3 Incorporating Asymptotic Freedom 

The goal of AdS / QCD is to reproduce a model that resembles QCD as close as possible, including 
QCD asymptotic freedom behaviour. Following Csaki et. al. formulation [42], we define our action 

^ ^ d^x^(-n-V{<p) + lg'^''dMCpdNcp). (3.1) 



with the 5-dimensional Newton constant and V{(p) chosen to reproduce asymptotic freedom. Note 
that (p is dimensionless here. 

We assume QCD gauge coupling is given by e^'?'(^) just like in string theory. For the coupling to 
run logarithmically we need a solution of the form: 

logf 

where as usual we have identified the AdS coordinate z with the energy scale and defined zg = ^qqd- 
We do not fix the value of a priori. It is easier to work in coordinate y given by exp (y/R) = z/R 
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in which our solution will take the form: e^'?'(y) = -g^t^- Inverting this we obtain the function for the 
dilaton 

and its derivative 

First we write our metric in the y coordinate: 

ds^ = e~^^^y\rif,ydxf'dx'' - dy^). (3.5) 
Subtituting the metric into our action we obtain: 

Sg^ = 2^ / ^'^V^ (-20^" + - ^(^^ - ■ (3.6) 

Taking the variance of this action with respect to A and (p and setting it to zero: 

SS 

= ^ 48A'^ - 24A" + W{(p) + l{(p'f = 0, 

f^ = 0^,"-4fA' = ^^ ,3.7, 

Next we use superpotential method by defining a function W(^) such that: 

A'{y) = W{cp{y)), 

cp'iy) = 6^. (3.8) 
Solving for W{(p) with (p' as defined in Eq. 3.4, we obtain: 

Wi^P) = ^e^^ + Wo- (3.9) 



It is very easy to show that the potential is 



y(^) = 18(-) -12Wl 



Now we can solve for the warp factor: 



(3.10) 



A(y) = Ao + Woy+^log^, (3.11) 

or in z coordinate: 

A(z) = Ao + WoRlog|-^loglog^. (3.12) 

For e^^'^^^i = e~^^° I ) (log y) 31^^ to resemble AdS metric with log Zq/z factor, we need Aq = 
and Wo = 1/R. This will give us 
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and 



For the case b = ± \/ 1, we got: 



and 



^ = T#og^ = TJfloglogf. (3.17) 



And finally the metric is 



1 

57 



2/1 \ 

log ^) ^ dxf'dx>,„, - . (3.18) 

In order to write down the action of the dilaton couple into matter in Einstein frame, we com- 
pared the action in string and Einstein frame of a noncritical string theory in 5 dimensions [42]. In 
string frame the action has the form [45, 46] 



2kI 



(3.19) 



S = ^ y d'x^ ^-n + -dM^d""'^ - V{^) + e-3'^C„j , (3.20) 
where the metric of the string frame and the metric of the Einstein frame are related as g^'f^ = 

4 ^ 

Smn- We can rescale the dilaton field to obtain a canonical kinetic term, by using <t> = \'3/8(p. 



The action in Einstein frame is: 



4 Vector Meson 

We have chosen the background dilaton profile such that the model produce asjonptotic freedom. 
Here we will calculate the resulting vector meson spectrum on this background. Using the action in 
equation (3.20) and the rescaled dilaton field <t> = ^/3/8(p, the vector part of the action can be written 
as 

Sv = l d^x^gei^^^ 1^^^ s'^'s'^'Vlp, (4.1) 

with the following equation of motion for the bulk-to-boundary propagator of the transverse vector 
field 

a^V(,,z) - f 1 + ~^f\!^ \ + ^^V(,,z) = 0, (4.2) 

where q denotes 4-momentum and h as defined in (3.2). 
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In order for the action evaluated on the solution of the equation (4.2) to be finite when b > 
(2/7) V2/3, the bulk-to-boundary propagator or its derivative must vanish sufficiently fast at the IR 
boundary. Multiplication of both should vanish sufficiently fast such that the IR surface term of the 
action vanish. 

Evaluating the action on the solution, we obtain 



s^ = /d^x^yO(x)prK°W^(inf)"^ 



V{q,z)d,V{q,z) 



(4.3) 



z=e->0 



where Pj^ — {rj^'" — qt'q^ /q^) and Vj^(x) is the UV boundary value of the vector field which corre- 
sponds to the source of vector current operator. In order to obtain 2-point function for the vector 
current, one has to take functional derivative with respect to the source Vjj twice on the 5D action. 
One obtains from the AdS / CFT correspondence 



Tr{x)f'{0) 



(4.4) 



where we have set the AdS radius to unity {R = 1). The bulk-to-boimdary propagator can be nor- 
malized to V(<j, e) = 1 at the UV boimdary. 

The vector meson mass spectrum can be obtained by solving the eigenvalue equation 



7 

"6F 



1 

W 



m 



Infi 



+ ^i/^„(z) = 0. 



(4.5) 
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where m„ is the mass of the n-th Kaluza Klein modes. The wave function is normalized as follows 



dzifln^ 

z V z 



7_ 3 

'66 V 2 



lpniz)^pm{z) = &mn- 



(4.6) 



The wave functions and their derivatives vanish at both the UV and the IR boundaries, except for the 
massless mode where the wave function approach a non-zero constant at IR boundary. The vanishing 
values of the wave functions at the boundaries is crucial for their normalizability. 
Setting mg = in equation (4.5), one obtains the following solution 



i/'o = NoT 1 




(4.7) 



which is the upper incomplete gamma function. This function is vanishing at UV boundary and 
approach Nor(l — 1/ (3^^) + 7/ (6^;) V3/2) at IR boundary. The slope vanishes at the IR boundary 
for h > (2/7) V2/3. The weighting function of the normalization integral is infinite at both UV and 
IR boundary. However, the vanishing value of ipo near UV overcome this to make the integrand in 
(4.6) goes to zero at the UV boundary. This is not the case at the IR boundary as the wave fimction 
goes to a non-zero constant when b > (2/7)-\/2/3. Near the IR boundary, one can pull out xpQ out of 
the integral. This leave us with the following integral 



dz- 

z 



In 



(4.8) 
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Figure 1. Left: Plot of vector meson's decay constant f j as a function of b. Right: Plot of Zq as a function of b. 
In both plot we fix the lightest vector meson's mass to 775.5 MeV 



which is infinite for < b < (7/6) ^/3/2, with ?/ close to zq. In the next section we will show that with 
b > the quark mass runs properly, that is, in the high energy limit as the energy scale increases, the 
quark mass decreases as in QCD. 



In the case of = \/2/3, the first four mass eigenvalues of equation (4.5) ( in unit of 1 / zg) are 
3.83, 6.52, 9.14, 11.74. We set the value of zq to fit the mass of the first mode to the mass of rho meson. 
For nip = 775.5 MeV, we obtain 1/zq = 202.6 MeV. The resulting mass of the first few excited states 
are 1321 MeV, 1852 MeV, 2379 MeV. This mass spectrum grows like mfj ^ n^, similar to the hard- 
wall model. As comparisons, in the original hard-wall model, the corresponding mass of the excited 
states are 1780 MeV, 2789 MeV, 3802 MeV. For the soft- wall model with quadratic dilaton profile, the 
corresponding mass are 1097 MeV, 1343 MeV, 1551 MeV 

The bulk-to-boundary propagator can be written as a sum over normalizable modes. 



This expansion can be substituted into (4.4) from which one can identify the decay constant of the 
vector meson 



where = In [13]. We obtain Fi = (222 MeV)'^. As comparisons, the hard-wall model and the soft- 
wall model with quadratic dilaton profile yield fi = (329 MeV)^ and fi = (260 MeV)^, respectively. 

We can vary the value of b while fixing the mass of the lightest vector meson to 775.5 MeV The 
results for the decay constant and the value of Zq^ are shown in Figure 1. 

5 Running quark mass and QCD condensate 

The action for dilaton coupled to matter in Einstein frame is: 




(4.9) 




(4.10) 




(5.1) 



with D^X = d^X - iAfX + iXA^ and F^^ = d^A^ 
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Figure 2. (Color online). Plot of running of the quark mass as a function of energy scale for different r. The blue 

( — ) is the running of quark mass from QCD. The green ( ) is for r = —0.29 {b = 0.52), the orange ( ) 

is for r = -0.53 (b = 0.60), and the red ( ) is for r = -0.75 {b = VV3). 



The equation of motion for the background field Xq in the u = (yo — y)/R coordinate: 

du (/"(M)''a.Xo) +3e4"(M)'-Xo = 0, 



(5.2) 



where r = 3I7 ~ ^ y § ■ The solutions of the equation of motion (5.2) from [47, 48] in terms of Kummer 
functions are: 



Xo(m) = Cie^3"M^"''M 



l + -,2-r,2M 



l + -,2-r,2u 
2 



(5.3) 



with Ci and C2 are integration constants. Contrary to the soft-wall model, we can keep both terms 
since the action evaluated on the above solution remains finite. 

In order to understand how the scalar field behaves near the UV boimdary, we write the solution 
in the following form, using coordinate transformation u = — ln(z/zo) and Kummer's transforma- 
tion [47, 48] 



Xo(z) = Ci(ln^°)i-'-LZ 



3'' ^ r,, z 
— ,2-/-,21n — 
2 zq 



■C2(ln 



r z 
-,2-r,-21n- 

2 zoj Vzo 



The constants are related as follows 



rf2 - r) 

Ci = -Ci ^,^ , _^/^^ exp{-i3nr/2), 



C2 = C2 - C; 



r(l + r/2) 
r(2-r) 



r(l-3r/2) 



exp(— z'7rr/2) . 



(5.4) 

(5.5) 
(5.6) 
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Figure 3. (Color online). Plot of running of the quark condensate as a function of energy scale for different r. 
The blue ( — ) is the running of quark mass from QCD. The green ( ) is for r = —0.29 {b = 0.52), the orange 
( ) is for r = -0.53 {b = 0.60) and , the red ( ) is for r = -0.75 {b = vWS). 



The Kummer's functions in equation (5.4) can be written as an infinite series 



+C2(2)-(l+'-/2) ± 



n=0 



2"m! 



In 21 



'3r^ 



' l)n \T)n 

2"n\ 



In 21 



(5.7) 



where {a)n = («)(« + 1) ... (a + ji — 1), and (fl)o = 1. In the limit z ^> both term in above equation 
vanish, however, the second term vanish faster than the first term because of the factor. Hence, in 
this limit we can drop the second term and identify the remaining term as quark mass multiplying 
as in (2.8). 

We will determine the constants of integration using the value of niq and ct" at 1/z = 2GeV 
from [49] as input parameter. Because the imaginary parts of equation (5.4) cancel out as a conse- 
quence of equation (5.3) which is real, we take the real part of Xo(z) and associate them to the quark 
mass nici and QCD condensate a: 



mcf{z) = Re 
(7(z) = Re 



Q(ln^)i-'-LI 
z 



C2(ln 



^°)i-LI 



l--,2-r,21n- 

2 zo 

r z 
1 + -,2-r,-2hl — 
2 Zo 



47r 
n 



(5.8) 
(5.9) 
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By associating zq with I/AqcD/ equation (5.8) and (5.9) show how the running of quark mass 
and QCD condensate as a function of energy scale 1 / z = f/. Figure (2) shows the running of the 
quark mass from equation (5.8) compared with the running of quark mass from perturbative QCD 
up to two loops calculation [50, 51], 



2j6i7olnL + l 



8_7i 



-^ + 0{L-^ln^L) 



(5.10) 



where = 2, = ^ — ^n^, /5o = 11 — /5i = 51 — ^"g, and L = ln(;;^/A^). Here, fhq is the 
renormalization group invariant mass. 

We vary the value of r while fixing the mass of the lightest vector meson nip = 757 MeV. The 
best fit as shown in figure (2) is for r = —0.53 which corresponds to I/zq = 226 MeV and y/T^ = 197 
MeV. 

In figure (3), we fit equation (5.9) to the one loop QCD calculation for the running of quark 
condensate in the chiral limit [52, 53]: 



A2 



(5.11) 



with 7m = 12/(33 — 2^,,). The best fit is for r = —0.29 which corresponds to I/zq = 259 MeV and 



/Fp = 183 MeV. 



6 Conclusion 



In this paper we have calculated the running of the quark mass and the quark condensate using a 
dynamical soft-wall model. In order to remove the massless vector meson and to make the quark 
mass rims properly we have determined the allowed value of parameter hioheO <b < (7/6)-\/3/2. 
We also find that, contrary to the soft-wall model with quadratic dilaton profile, the model allows 
for a non-vanishing quark condensate in the chiral limit. We obtain the vector meson mass spectrum 
similar to that of the hard-wall model. The best fit of this model to the running of the quark mass 
yield r = —0.53, which corresponds ioh = 0.60 and the best fit to the running of the quark condensate 
yield r = —0.29, which corresponds ioh — 0.52. 
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